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Recently superconductivity at the interface between the insulators LaAlOs and SrTiOa has been 
tuned with the electric field eff'ect to an unprecedented range of transition temperatures. Here we 
perform a detailed finite size scaling analysis to explore the compatibility of the phase transition line 
with Berezinskii-Kosterlitz-Thouless (BKT) behavior and a 2D-quantum phase(QP)-transition. In 
an intermediate regime, limited by a gate voltage dependent limiting length, we uncover remarkable 
consistency with a BKT-critical line ending at a metallic quantum critical point, separating a weakly 
localized insulator from the superconducting phase. Our estimates for the critical exponents of the 
2D-QP-transition, z ~ 1 and 77 ~ 2/3, suggest that it belongs to the 3D-xy universality class. 
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I. INTRODUCTION 

At the interface between oxides, electronic properties 
have been generated, different from those of the con- 
stituent materials .-ii^iS, In particular, the interface be- 
tween LaAlOs and SrTiOs, two excellent band insula- 
tors, found to be conducting in 2004^ attracted a lot 
of attentioii^'^'^'i'^'^. Recently, different ground states, 
superconducting and ferromagnetic, have been reported 
for this fascinating systemi^ In a recent repori^iS, it was 
shown that the electric field effect can be used to map 
the phase diagram of this interface system revealing, de- 
pending on the doping level, a superconducting and non- 
superconducting ground state and evidence for a quan- 
tum phase transition. 

Continuous quantum phase transitions are transi- 
tions at absolute zero in which the ground state of 
a system is changed by varying a parameter of the 
Hamiltonian.— li^ii^ The transitions between supercon- 
ducting and insulating behavior in two-dimensional sys- 
tems tuned by disorder, film thickness, magnetic field or 
with the electrostatic field effect are believed to be such 
transitions^^'ii^'i^^^^ 

Here we present a detailed finite size scaling analysis 
of the temperature and gate voltage dependent resistivity 
data of Caviglia et alJ^ to explore in the LaAlOs/SrTiOs 
system the nature of the phase transition line and of its 
endpoint, separating the superconducting from the in- 
sulating ground state. For this purpose we explore the 
compatibility of the normal state to superconductor tran- 
sition with Berezinskii-Kosterlitz-Thouless (BKT) criti- 
cal behavior Our analysis of the temperature depen- 
dence of the sheet resistance at various fixed gate voltages 
uncovers a rounded BKT-transition. The rounding turns 
out to be fully consistent with a standard finite size effect 
whereupon the correlation length is prevented to grow 
beyond a limiting length L. Indeed, a finite extent of the 
homogeneous domains will prevent the correlation or lo- 
calization length to grow beyond a limiting length L and, 
as a result, a finite size effect occurs. Because the cor- 



relation length does not exhibit the usual and relatively 
slow algebraic divergence as Tc is approached, the BKT- 
transition is particularly susceptible to such finite size 
effects. Nevertheless, for sufficiently large L the critical 
regime can be attained and a finite size scaling analysis 
provides good approximations for the limit of fundamen- 
tal interest, L oori^'^^^ 

As will be shown below, our finite size scaling analysis 
uncovers close to the QP-transition a gate voltage de- 
pendent limiting length. According to this electrostatic 
tuning does not change the carrier density only but the 
inhomogeneity landscape as well. The finite size scaling 
analysis also allows us to determine the gate voltage de- 
pendence of the BKT-transition temperature Tc, of the 
associated fictitious infinite system. This critical line, Tc 
versus gate voltage, ends at a quantum critical point at 
the gate voltage Vgc- Here the sheet conductivity tends to 
(To {T ^ 0, Vgc) ~ 2.52 • 10"'' (fi-i) which is comparable 
to the quantum unit of conductivity Ae^ /h ~ 1.55 • 10~^ 
(il^^) for electron pairs, emphasizing the importance of 
quantum effects. Its limiting temperature dependence 
points to Fermi liquid behavior at quantum criticality. 
The estimates for the critical exponents of the 2D-QP- 
transition, z ~ 1 and F ~ 2/3, suggest that it belongs to 
the 3D-xy universality class. In the normal state we ob- 
serve non-Drude behavior, consistent with the evidence 
for weak localization. To identify the nature of the in- 
sulating phase from the temperature dependence of the 
resistance, we perform a finite size scaling analysis, re- 
vealing that the growth of the diverging length associated 
with weak localization is limited and gate voltage depen- 
dent as well. Nevertheless, we observe in both, the tem- 
perature and magnetic field dependence of the resistance, 
the characteristic weak localization behavior, pointing to 
a renormalized Fermi liquid. In addition we explore the 
Tc dependence of the vortex core radius and the vortex 
energy. These properties appear to be basic ingredients 
to understand the variation of Tc- In the superconducting 
phase we observe consistency with the standard quantum 
scaling form for the resistance, while in the weakly local- 
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izcd phase it appears to fail. In contrast to the quantum 
scahng approach we obtain the scahng function in the 
superconducting phase exphcitly. It is controlled by the 
BKT-phase transition line and the vortex energy. 

In Section II we sketch the theoretical background and 
present the detailed analysis of the resistivity data of 
Caviglia et alr^ We close with a brief summary and some 
discussion. 



II. THEORETICAL BACKGROUND AND DATA 
ANALYSIS 

A. BKT-TRANSITION 

To explore the compatibility with BKT critical behav- 
ior we invoke the characteristic temperature dependence 
of the correlation length above Tcr^ 

e (T) = Co exp f 27r/ (bt'/A ) , t = \T/T, - 1| , (1) 



where is the classical vortex core radius and b is re- 
lated to the energy needed to create a vortex i^^i^^i^^i^^ 
Note that b also enters the temperature dependence of 
the magnetic penetration depth A below the universal 
Nelson-Kostcrlitz jumi>2i: 



(T,)/A2(T)= (l + 5|t|^/V4) 



(2) 



Moreover, b is related to the vortex energy Ec in terms 

g £26.28 



b = f {E,/ (kBT,)) ^ 



(3) 



Invoking dynamic scaling the resistance R scales in 13 = 2 



RocC 



(4) 



where Zd is the dynamic critical exponent of the classical 
dynamics. Zd is usually not questioned to be anything 
but the value that describes simple diffusion: Zd ~ 2^ 
Combining these scaling forms we obtain 



with 



AirT^^yb, Ro K 1/^0- 



(6) 



Accordingly the compatibility of experimental resistivity 
data with the characteristic BKT-behavior can be ex- 
plored in terms of 



{d\nR/dT) 



-2/3 _ 



(7) 



Because the correlation length does not exhibit the usual 
and relatively slow algebraic divergence as T^, is ap- 
proached (Eq. ([T])) the BKT-transition is particularly 



susceptible to the finite size effect. It prevents the cor- 
relation length to grow beyond a limiting lateral length 
L and leads to a rounded BKT-transition. Nevertheless, 
for sufficiently large L the critical regime can be attained 
and a finite size scaling analysis allows good approxima- 
tions to be obtained for the limit L — > including 
estimates for Tc, bn, Rq, and their gate voltage depen- 
dence. In the present case potential candidates for a lim- 
iting length include the finite extent of the homogenous 
regions and the failure to cool the electron gas down to 
the lowest temperatures. In the latter case L is given 
by the value of the correlation length at the temperature 
where the failure of cooling sets in. In any case finite size 
scaling predicts that R (T, L) adopts the form 



R (T, L) 
R {T, oo) 



e(T,(x)) 

R{T,L) 



where 



exp 



Ro 



b„ IT - 



exp fefl |T - Tc 



-1/2 



1-1/2 



(8) 



(9) 



g(x) is the finite size scaling function. If ^ (T, oo) < L 
critical behavior can be observed as long as g(x) ~ 1, 
while for ^ (T, oo) > L the scaling function approaches 

g{x) oc X so R{T)exp(bR\T -Tc\^^^^^ /Ro tends to 

ig/Ro)cxp(bR \T - Tel"'/') with g ^ 1/L\ 

We are now prepared to explore the evidence for 
BKT-behavior. In Fig. [T] we show {dlnR/dTy'^^^ 
vs. T for = 40 V. In spite of the rounded tran- 
sition there is an intermediate regime revealing the 
characteristic BKT-behavior ([7]), allowing us to esti- 
mate i?o, bfi and T^,. As can be seen in the in- 
set of[Il depicting R{T) cxp(^5h |T - TJ"^/^^ /i?o vs. 

(l/i?o)exp|^6_R |r — Tcl"^''^^ , the rounding of the tran- 
sition is remarkably consistent with a standard finite size 
effect. The horizontal line corresponds to ^ < L where 
critical behavior can be observed as long as g{x) ~ 1, 
while the dashed one characterizes the rounded regime 
where ^ > L. Here the scaling function approaches 
^(a::) oc x and R{T,L) tends to <? oc L~^. Independent 
evidence for BKT-behavior was also established in ear- 
lier work in terms of the current-voltage characteristics!^ 
Applying this approach to the R (T) data for each gate 
voltage Vg we obtain good approximations for the val- 
ues of Tc{Vg), bR{Vg), and Ro{Vg), in the absence of 
a finite size effect. The resulting BKT-transition line is 
depicted in Fig. displayed as Tc vs. i?n (T*), the 
normal state resistance at T* ~ 0.4K. We observe that 
it ends around R^jc (T*) — 4.28 kfi where the system 
is expected to undergo a 2D-QP-transition because Tc 
vanishes. With reduced R\j the transition temperature 
increases and reaches its maximum value, Tcm — 0.31 K, 
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FIG. 1: (dlni?/rfr)~^/^ vs. T for Vg = 40 V where R = 
3/5Rn. The sohd hne is [din R/dT)~^/^ = 6.5 (T - 
yielding the estimates Tc = 0.27 K and (2/6^)^/^ = 
6.5; the inset shows R{T) exp(bR \T - Tc^'^^^^ / Rq vs. 

{l/Rt,)exp(bR \T - Tcr^/^) with i?o = 1-67 k Q. The upper 

branch corresponds to T > Tc and the lower one to T < Tc. 
The solid line is R {T, L) ~ R {T, oo) and the dashed one 

i? exp^&Hlr-Tcr'''') /i?o = ig/Ro)exp(bR\T -Tc]-''^^^ 

with p ~ 501 (X 

around i?n (T*) ~ 1.35 kfi. With further reduced resis- 
tance Tc decreases. We also included the gate voltage 
dependence of the normal state resistance since correc- 
tions to Drudc behavior (a (x n) have been discussed in 
the literature for systems exhibiting weak localization as 
will be demonstrated below.— 




R (T) (n) 

FIG. 2: Tc vs. Rn (T*) (•) and Vg vs. Rn (T*) (*) at T* = 0.4 
K. The error bars indicate the uncertainty in the finite size 
estimates of Tc. The solid line is Tc = 1.17 x 10"*Ai?(r*) 
and the dashed one Vg = Vgc + 1.39 x 10"^ A7?^/2(T*) with 
Aii(T*) = {Ruc (T*) - Rn (T*)), Roc {T") = 4.28 kfi and 
Vgc = -140 V. 

According to the scaling theory of quantum criti- 
cal phenomena one expects that close to the 2D-QP- 



transition Tc scales a.fM^ 

Tc oc <5^^ (10) 

where 6 is the appropriate scaling argument, measuring 
the relative distance from criticality. V denotes the crit- 
ical exponent of the zero temperature correlation length 
^ (T = 0) oc (S"" and z the dynamic critical exponent. 
From Fig. [2] it is seen that the experimental data points 
to the relationship 

T,cx Ai?n(T*)« (11) 

close to quantum criticality, where Ai?n (T*) = 
RadT*) - Rn (T*). In this context it is important to 
emphasize that Tc cx Ai?n (T*) turns out to be nearly 
independent of the choice of T* around T* « 0.4 K. So 
the normal state sheet resistance R\j (T*) is an appropri- 
ate scaling variable in terms of Ai?n (T*). In this case 
zV = I, while if (5 = AVg, zV — 2/3. Since the measured 
modulation of the gate voltage induced charge density 
An2D scales in the regime of interest asi^ 

AVg OC An2D cx T^/^, (12) 

so zV = 2/3 if AVg or An2D are taken as scaling argu- 
ment S. On the other hand it is known that S oc An2D 
holds if {2 + z)V > 2^ To check this inequality, given 
zV, we need an estimate of z. For this purpose we invoke 
the relation Rq — i?oc oc ^o"^ (Eq- ®) and note that the 
critical amplitude of the finite temperature correlation 
length ^0 and its zero temperature counterpart should 
scale as c< ^(T = 0) oc 6~'^, so that the scaling relation 

Ra - Rac oc Co^' oc (T = 0) oc (5^^ oc (13) 

holds. Fig. [3] depicts the Tc dependence of the vor- 
tex core radius oc ^ (T = 0) oc (i?oc — ^o) and 
b, which is related to the vortex energy Ec- Approach- 
ing the 2D-QP-transition we observe that the data point 
to ^ (T = 0) oc 1/Tc, yielding for z the estimate z ~ 1 
so that P 2/3 with zD ~ 2/3. As these exponents 
satisfy the inequality {2 + z)V > 3^ we identified the 
correct scaling argument. S oc An2D oc AVg. The 2D- 
QP-transition is then characterized by the scaling rela- 
tions 

Tc oc oc ARa (T*) oc AF//^ oc Anl^j^ oc (14) 

where ARu (T*) oc An^^^ reveals non-Drude behavior 
in the normal state. The product zV ~ 2/3 agrees 
with that found in the electric field effect tuned 2D- 
QP-transition in amorphous ultrathin bismuth filmsi^ 
and the magnctic-field-induced 2D-QP transition in 
Nbo.15Sio.85 filmsi^ On the contrary it differs from the 
value zF ~ 1 that has been found in thin NdBa2Cu307 
films using the electric-field-effect modulation of the tran- 
sition temperatureii^ In any case our estimates, z ~ 1 
and P ~ 2/3 point to a 2D-QP-transition which belongs 
to the 3D-xy universality classii^ 
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FIG. 3: Vortex radius « {R'oc - Ro)~^^^ (•) and b (*) vs. 

where R = 3/5Rn . The soUd hne is Co <x {Roc - Ro)~^^^ = 
8 X 10"^/Tc with Roa = 2.7 kfJ. 
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FIG. 4: Gate vohage dependence of the limiting length L in 
terms of L cx g~^^^ vs. Vg. The inset shows Tc vs. Vg. The 
solid line is Tc = 8.9 x 10~^{Vg - Vg^f^^ (K) indicating the 
leading quantum critical behavior (|14|l with zu ^ 2/3. 



Fig. [3] also depicts the Tc dependence of b, which is 
related to the vortex energy. Since b tends to a constant 
in the hmit 0, Eq. (O imphes db/dT^ = and 

therewith 

E,{T,) oc keT,, (15) 
while the core radius diverges as 

?o K i/r„ (16) 

in analogy to the behavior of superfluid '*He films where 
Tc was tuned by varying the film thickness.— A lin- 
ear relationship between the vortex core energy and 
Tc was also predicted for heavily underdoped cuprate 
superconductorsj^ Furthermore, an increase of the vor- 
tex core radius with reduced Tc was also observed in 
underdoped YBa2Cu30j;^ and La2-a;Sra;Cu04i2i The 
2D-QP-transition is then also characterized by vortices 
having an infinite radius and vanishing core energy. As 
Tc increases from the 2D-QP transition, the core radius 
shrinks, while the vortex energy increases. We also ob- 
serve that the rise of Tc is limited by a critical value of 
the core radius and that the maximum Tc {Tcm — 0.31 K) 
is distinguished by an infinite slope of both, the vortex 
radius and b. Finally, after passing Tcm the vortex core 
radius continues to decrease with reduced Tc while b 
increases further. 

Next we explore the gate voltage dependence of the 
limiting length. Indeed, its presence or absence al- 
lows us to discriminate between an intrinsic or ex- 
trinsic limiting length. For this purpose we per- 
formed the finite scaling analysis outlined in Fig. 1 
for various gate voltages. In the finite size domi- 
nated regime, ^ > L, the finite size scaling form ([8]) 

reduces to R{T)exp(bB.\T - Tc\~^^^^ / Rq = g{x) oc 

{g/Ro)cxp(^bR \T — Tcl""^^^^ with g oc 1/L^, so g probes, 
if there is any, the gate voltage dependence of L. In 



Fig. 4 we summarized the resulting gate voltage depen- 
dence of L oc g~^^^. It is seen that the limiting length is 
nearly gate voltage independent down to = V. This 
points to the presence of inhomogeneities preventing the 
correlation length to grow beyond the lateral extent of 
the homogeneous domains. On the contrary, for nega- 
tive gate voltages L decreases by approaching the QP- 
transition as Tc does. The resulting broadening of the 
BKT-transition with reduced Vg and Tc is apparent in 
the temperature dependence of the sheet resistances^ A 
potential candidate for a gate voltage dependent limiting 
length is the failure of cooling at very low temperatures 
In this case the correlation length cannot grow beyond its 
value at the temperature Tf where the failure of cooling 
sets in. Invoking Eq. ([T|) in the limit Tc ^ we obtain 

Lf = ^ocxp^27r/(6T^^^^^ oc 1/Tc, because c>c 1/Tc and 

b remains finite in the limit Tc — > (see Fig. [3]). Con- 
trariwise we observe in Fig. 2] that L decreases with Tc. 
According to this electrostatic tuning docs not change the 
carrier density only but the inhomogeneity landscape as 
well. 

In any case, the agreement with BKT-behavior, lim- 
ited by a standard finite size effect, allows us to dis- 
criminate the rounded transition from other scenarios, 
including strong disorder which destroys the BKT- be- 
havior. It also provides the basis to estimate Tc, bu 

— 1/9 

and i?o, and with that b and oc i?g with reason- 
able accuracy. The resulting BKT-transition line ends 
at Vgc — —140 V, where Tc vanishes and the system 
undergoes a 2D-QP transition. After passing this tran- 
sition Tc increases with reduced negative gate voltage, 
reaches its maximum, Tc„i — 0.31 K, around Vg ~ 110 V 
and decreases with further increase of the positive gate 
voltage. Remarkably enough, this uncovers a close anal- 
ogy to the doping dependence of Tc in a variety of bulk 
cuprate superconductorS j^^i"^^ where after passing the so 
called underdoped limit Tc reaches its maximum with in- 



5 



creasing dopant concentration. With further increase of 
the dopant concentration decreases and finally van- 
ishes in the overdoped limit. This phase transition line 
is thought to be a generic property of bulk cuprate su- 
perconductors. There is, however, an essential difference. 
Cuprates are bulk superconductors and the approach to 
the underdopcd limit, where the QP transition occurs, is 
associated with a 3D to 2D crossover fi^*^ while in the 
present case the system is and remains 2D, as the consis- 
tency with BKT critical behavior reveals. Furthermore, 
a superconducting dome (in the T versus doping phase 
diagram) was also observed in bulk doped SrTi03 that 
is close to the system under studyj^i^ 

B. INSULATING PHASE 

Supposing that the insulating phase is a weakly local- 
ized Fermi liquid the sheet conductivity should scale as^ 

an(r) =aao+dln(r), (17) 

where d = e^/ (tt/i) ~ 1.23 x 10"^ is generically 

attributed to electron-electron interaction;^ while crno 
is expected to depend on the gate voltage. In Fig. [5^ 
we depicted ctq — ctoo vs. T for various gate voltages 
Vg by adjusting to achieve a data collapse at suf- 
ficiently high temperatures. The resulting gate voltage 
dependence of ctdo, consistent with 

^□o(K,) = an, -5.9x10-^ IK,- t/<,,|'/^(f}-i), 
(JOs = 2.52 X (18) 

is shown in Fig. [5)d. An important feature of the data 
is the consistency with a weakly localized Fermi liquid 
because the coefficient d is close to d = e^/ (tt/i). In 
any case, more extended evidence for weak localization 
emerges from the magnctoconductivity presented below 
(Fig.®. 

A very distinct temperature dependence of the conduc- 
tivity occurs at quantum criticality, Vg = —140 V~ Vgc- 
Indeed, the dashed line in Fig. [5^ and the solid one in 
Fig. [6] indicate that in the limit T — *■ the system tends 
towards a critical value. According to the plots shown in 
Fig. [5] the limiting behavior is well described by 

an (T, Vgc) = au, {Vgc) - 9.782 x IQ-^T^. (19) 

Note that our estimate <Tn (T = 0, Vgc) = o-qo {Vgc) ~ 
aus = 2.52 X 10~^J7~^ is comparable to the quantum 
unit of conductivity Ae'^/h ~ 1.55 x 10~^ rt~^ for elec- 
tron pairs, emphasizing the importance of quantum ef- 
fects. The dependence points to Fermi liquid behavior 
in the regime ksT « hiojj, Ep where electron-electron 
scattering dominates, lod is the Debye frequency and 
Ep denotes the Fermi energy. At higher temperature 
we observe a crossover to a linear T -dependent con- 
ductivity marked by the dash-dot line. Recent theories 
on the conductivity of 2D Fermi liquids predict such a 
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FIG. 5: a) crn — a^o 'vs. T for various gate voltages Vg. The 
solid line is an — crno ~ d ln(r) (fi^^) with d ~ 1.23 x 
10^'' and (Tno (Vg) taken from Fig. [5j3. The dot at the 
origin marks the quantum critical point and the dashed line 
is Eq. (|19p indicating the low-temperature behavior of the 
sheet conductivity at the quantum critical point, b) ctqo vs. 
Vg. The solid line is Eq. with Vgc = -140 V. 

linear r-dependence4^ From Eqs. pT)) and (HH]), de- 
scribing the data in the weakly localized regime rather 
well, it also follows that the normal state conductiv- 
ity at T* = 0.4 K scales as aoc {T*) - an {T*) cx 
10~^ \Vg — Vgcf^^. Together with the empirical scahng 
relation pT)) . \Vg — Vgd cx An2c it points to non-Drude 
behavior in the normal state. 

Considering the temperature dependence of the sheet 
conductivity below quantum criticality {Vg = —140 V~ 
Vgc), Fig. [S^ reveals at sufficiently high temperature re- 
markable agreement with the ln(r) behavior, character- 
istic for weak localization. On the contrary, in the low 
temperature regime and even rather deep in the insulat- 
ing phase {Vg = —300 V), systematic deviations occur in 
terms of saturation and an upturn as quantum criticality 
is approached {Vgc = —140 V). Because the conductivity 
of a weakly localized insulator is not expected to saturate 
in the zero temperature limit^i^ this behavior appears 
to be a finite size effect, preventing the diverging length 
associated with localization^, ^,oc oc (i|ln(T)|, to grow 
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T(K) 

FIG. 6: aa {Vgc,T) vs. T and in the inset vs. at Vgc = 
— 140 V. The solid line, indicating consistency with is Eq. 
(fT9)) . while the dashed line is an {Vgc,T) = 2.275 x lO""* + 
1.54 X 10"^r (r2"^). The the dash-dot one in the inset is 
again Eq. (fT9)l . 
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FIG. 7: (an (T) -ano)/(rfIn(r)) vs. 1/ {d\ln{T)\) at Vg = 
-220, -240 and -300 V with d = 1.23 x 10"^ The lines 

correspond to Eq. (|2ip providing a measure for L in terms of 
\gL\ oc L. 



beyond L, the limiting length already identified in the 
context of the rounded BKT-transition (Fig. 01). In the 
present case finite size scaling predicts that a\j (T) should 
scale as 



□ c 



dln(r) 



5(y), y-L/6ocCxL/(d|ln(r)|). 



(20) 

g{y) is the finite size scaling function which tends to 1 
for y < 1. In this case the approach to the insulating 
ground state can be seen, while for y > 1 the crossover 
to g{y) — > y sets in and (Tn(T) approaches the finite size 
dominated regime, where 



o'njT) - g-pe 
filn (T) 



= gL/{d\\n{T)\), gL^L (21) 



A glance at Fig. El depicting {aa (T) - aao) / (dln(T)) 
vs. l/(dln(T)) at Vg = -220, -240 and -300 V, re- 
veals that the systematic deviations from the character- 
istic weak localization temperature dependence are fully 
consistent with a standard finite size effect. Accordingly, 
the saturation and upturns seen in Fig. [5] at low tem- 
peratures are attributable to a finite size effect, while 
in a homogeneous and infinite system the data should 
collapse on the solid line in Fig. [5^. An essential excep- 
tion is = —140 V. Here the interface approaches the 
metallic quantum critical point (see Fig. [5]) , metallic be- 
cause the sheet conductivity remains finite, approaching 
an (T = 0, Vgc) - <Jao {Vgc) = eras ^ 2.52 x lO"" 



(Eq. IHl)) in the limit T 0. 

Fig. 7 also reveals that the limiting length, L oc g^, de- 
pends in the insulating phase on the gate- voltage as well. 
The resulting dependence, L (Vg) oc g^ (Vg), is shown in 
Fig. [H In analogy to the limiting length associated with 
the BKT-transition (Fig. 01) it decreases by approaching 



quantum criticality at Vg 



-140 V. As a reduction of 



feature accounts for the saturation and upturns seen in 
Fig. [5^. Supposing that the limiting length is set by the 
failure of cooling below the temperature Tf then L is set 
by Ly = £,ioc{Tf) oc (i|ln(T/)| and with that indepen- 
dent of the gate voltage, in disagreement with Fig. [S] 
Accordingly, in analogy to the BKT-transition, the lim- 
iting length appears to be attributable to a electrostatic 
mediated change of the inhomogeneity landscape. 




L enhances deviations from the asymptotic behavior this 



150 



FIG. 8: —gL oc L vs. Vg derived from finite size scaling plots 
as shown in Fig. [7] with Eq. (PT|) . 

Direct experimental evidence for a limiting length 
emerges from the work of Hani et al^ A single electron 
transistor was used as a local electrostatic probe to study 
the underlying spatial structure of the metal-insulator 
transition in two dimensions. The measurements show 
that as the transition is approached from the metallic 
side, a new phase emerges that consists of weakly cou- 
pled fragments of the two-dimensional system. These 
fragments consist of localized charge that coexists with 
the surrounding metallic phase. As the density is low- 
ered into the insulating phase, the number of fragments 
increases on account of the disappearing metallic phase. 
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The measurements suggest that the metal-insulator tran- 
sition is a result of the microscopic restructuring that oc- 
curs in the system. On the other hand, we have seen that 
the limiting length associated with the resulting inhomo- 
geneities depends on the gate voltage (see Figs. |4] and 
El). 

Further evidence for a weakly localized insu- 
lating phase stems from the observed negative 
magnetoresistanceiiS An applied magnetic field leads 
to a new length given by the size of the first Landau 

1 /2 

orbit, or magnetic length, Lh — {^q/ {2ttH)) , which 
decreases with growing field strength. Once its size 
becomes comparable to the dephasing length Lth (dis- 
tance between inelastic collisions)^ weak localization 
is suppressed. In D=2 the following formula for the 
magnetoconductivity was obtainedi ^^i^^ 



f^n = ffno + c [V- (1/2 + l/x) + In (x)] , c : 



Th 



(22) 



where ip denotes the digamma function, a* is a constant 
of the order of unity,— and 



X 







In the limit a; >> 1 it reduces to 

,2 



a e 



nh 



-1.96 + In (x)] 



while in the limit x 



holds. Here 



dag 
d\n{H) 



a*1.24 X IQ-^n 



-5o-i 



(23) 



(24) 



(25) 



(26) 



applies. In Fig. [9] we compare the experimental data with 
the theoretical predictions. The data agrees reasonably 
well with the characteristic weak localization behavior 
(Eq. (|22)) ). while the asymptotic In(iJ) behavior (Eq. 
([M)) ) is not fully attained. The resulting estimates for 
daa/d\n{H) are close to e^/ (tt/i) ~ 1.24 x IQ-^rj^i and 
consistent with the zero field temperature dependence 
of the sheet conductivity, ag (T) = ago + d\n(T), with 
d ~ e^/ (tt/i) (see Fig. [S^). An analogous treatment 
of the magnetoresistance data of a non superconducting 
sample of Brinkman et alA yields 87riy^/$o — 2.862 
T-i and c = a*e'^/TTh = 4.8 x 10"^ so c adopts 

in 'superconducting' and 'non superconducting' samples 
substantially different values. In any case, our analysis 
of the magnetoconductivity uncovers a weakly localized 
insulating phase, consistent with the ln(T) temperature 
dependence of the zero field counterpart at sufficiently 
high temperatures, and non-Drudc behavior in the nor- 
mal state. 
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FIG. 9: Magnetoconductivity erg vs. H, applied perpendic- 
ular to the interface, at T = 0.03 K and Vg = —300 V and 
-340 V. The solid line is an = 4.51 x 10"^ + 1.2 ■ 10"^ In (H) 
kSl"\ the dashed one ctq = 2.6 x 10"^ -I- 1.1 x 10'^\n {H) 
kfl~^, the dotted and dash dot curves are Eq. p2|) with the 
aao, 87rL|,^/$o and c = Q*e^/(7r/i) values 0.0284 kr2"\ 2.906 
T-\ 1.46 x lO"'^ Q-^ for Vg = -340 V and 0.044 kQ'^ 



4.068 T-\ 1.46 X 10"^ fi"^ at Va 



-300 V. Note that 



87rL|^/$o = 2.906 T"^ corresponds to Lth ^ 1.55x 10"®cm 
whereupon Lh = Lth at H — 1.37 T. 



C. QUANTUM PHASE TRANSITION 

Traditionally the interpretation of experimental data 
taken close to the 2D-QP-transition was based on the 
quantum scaling relatio n^"^'^^i^° 



Rn(T,S) = BasG{x), x = cS/T 



l/z 



(27) 



G {x) is a scaling function of its argument and G (0) = 1, 
so at quantum criticality the system is metallic with sheet 
resistance Ros- The BKT-line is then fixed by Xc = 
cS/Tc^'"', whereby vanishes as ex oc ^"^ (T = 0) 
(Eq. P^ ). c is a nonuniversal parameter and S the ap- 
propriate scaling argument, measuring the relative dis- 
tance from criticality. This scaling form follows by not- 
ing that the divergence of ^ (T = 0) oc (5"" is at finite 
temperature cut off by a length Lt, which is determined 
by the temperature: Lt (x T~^^^. Thus G (x) is a finite 
size scaling function because x oc {Lt/£, {T = O))^^'' oc 

The data for Rn{T,5) plotted vs. 
should then collapse onto two branches joining at R[2s- 
The lower branch stems from the superconducting and 
the upper one from the insulating phase. To explore the 
consistency with the critical BKT-behavior we note that 
in the limit ^ the relation 



Rkt iT,Vg) 



exp 



iT-T4Vg)) 



1/2 



Gl'Y^ 



(28) 



should apply. Indeed, the BKT-scaling form of the re- 
sistance applies for any T ^ T^^ because the univer- 
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sal critical behavior close to Tc is entirely classical)^ 
On the contrary T^, bn and the critical amplitude Rq 
arc non-universal quantities which depend on the tun- 
ing parameter. Furthermore, they are renormalized by 
quantum fluctuations. In any case, the data plotted as 
Rkt (T, Vg) /Ro (Vg) vs. {Vg - Vgc) /T^^''' should col- 
lapse on a single curve and approach one close to quan- 
tum criticality. In Fig. [10] wc depicted this scaling plot, 
derived from Rq (Vg), Ir {Vg) and {Vg) for zV ^ 2/3. 
Apparently, the flow to the quantum critical point is well 
confirmed. Furthermore, noting that close to the QP- 



and m]) as 



transition 6^ (V^) oc Tc^^, because = A^nTc'^/b (Eq. 
dH])) and b ~ const, (see Fig. [3]), the vortex core energy 
scales as Ec oc ksTc (Eq. (fT5|) ). the scaling function 
adopts with zD 
the form 



1/2, 



2/3 and oc {Vg - Vgcf^^ (Eq. (HI 



G{x) ~ exp(-Sxi/V(l - fea;2/3)i/2)^ ^ 



shown by the solid line in Fig. (TUl 



Since Vg - Vgc 



(29) 
oc 



AttTc' /b is related to the vor- 



Tc^^'^ oc Tc^"^ and b^ 
tex energy in terms of b (Eq. ([3])), the scaling function 
is controlled by the BKT-line and the vortex core en- 
ergy, while the vortex core radius enters the prefactor 
via Ro {Vg) - Roc oc ^o"' oc (T = 0) oc (Eq. ^). 
Noting that a ~ Ana^^^/b and b = a, where a is given 
in terms of Tc = a{Vg - Vgc)^/^ (K) with a ~ 8.9 • lO^^ 
(see Fig. 15) and 6 ~ 50 (Fig. ^ we obtain a ~ 0.0237 
and 6 = a ~ 0.0089, in reasonable agreement with the 
fit parameters yielding the solid line in Fig. (TUl This 
uncovers the consistency and reliability of our estimates 
along the BKT-line. In this context it should be kept in 
mind that our analysis of the insulating state is limited 
by the finite size effect, preventing to approach the zero 
temperature regime. 




200 400 600 

>'=(V-V,)/T=")(VK-'") 

FIG. 10: RKT{Vg,T)/Ro{Vg) vs. {Vg - V,c) /T'^^ for var- 
ious Vg's. The solid line is Eq. ^ with a = 0.0248 and 
& = 0.0081 in K^^^VK 

On the contrary, in the insulating phase we observed 
that the sheet conductivity scales according to Eqs. (fT7|) 



= 1 - '-^^ \Vg - Vgc^ + — ln(T). 

(30) 

which is incompatible with the standard scaling form 
(|27p . Indeed, it involves two independent lengths. 
^id oc l/|ln(r)|, the diverging length associated with 

localization^ and ({T ^0) x \AVg\^^^'^ , the zero tem- 
perature correlation length (Eq. (|14p ). In this context it 
should be kept in mind that our analysis of the insulating 
state does not extend to zero temperature because Eq. 
([50)) applies at finite temperatures only. As T is reduced 
further the question of what happens in the insulating 
phase remains. 

To complete the BKT- and 2D-QP-transition scenario 
measurements of the magnetic penetration depth, A (T), 
would be required. At the BKT- transition Tc and A (T) 
are related by 



Ps {Tc) 



d$2 



167r3A2 {Tc) TT 



knTr. 



(31) 



while Ps {T) = above Tc. Ps is the 2D superfluid 
density and d is the thickness of the superconducting 
slieet^. The presence or absence of the resulting Nelson- 
Kosterlitz jump would then allow to discriminate experi- 
mentally between weak and strong disorder. In this con- 
text wc note that there is the Harris criterion,— which 
states that short-range correlated and uncorrelated dis- 
order is irrelevant at the unperturbed critical point, pro- 
vided that 2 — Diy < 0, where D is the dimensionality of 
the system and v the critical exponent of the finite tem- 
perature correlation length. With D = 2 and ly = oo, 
appropriate for the BKT-transition)^ any rounding of 
the jump should then be attributable to the finite size ef- 
fect stemming from the limiting length L. Furthermore, 
there is the quantum counterpart of the Nelson-Kostcrlitz 
relation, stating that 



d 



A2 (T = 0) 



16Tr^kBTcQ2 



(32) 



close to the 2D-QP transition.— li^i^ Q2 is a dimension- 
less critical amplitude bounded by^ 



< O2 < 1.11. 



(33) 



The lower bound corresponds to the BKT-line, 
d/A^ (T 0) ~ 1.03Tc with d, A in cm and Tc in K. Be- 
low this line the superfiuid order would become unstable 
to unbinding of vortices. The upper bound, corresponds 
to d/A^ (T = 0) ~ l.eiTc and the transition at T = be- 
longs to the BKT-universality class and consequently at 
Tc the superfluid density exhibits the universal disconti- 
nuity. Given our evidence for a (2-|-l)-xy QP transition, 
quantum fluctuations are present and expected to reduce 
Q2 from its maximum value, while the finite-temperature 
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transition remains again in the BKT universality classi^ 
Correspondingly, measurements of the temperature and 
gate voltage dependence of the supcrfluid density would 
be desirable to explore the observed BKT-behavior, weak 
localization and Fermi liquid features further. 

III. SUMMARY AND DISCUSSION 

In summary, we have shown that the electrostatically 
tuned phase transition line at the LaAlOs/SrTiOs in- 
terface, observed by Caviglia et al.,— is consistent with 
a BKT-line ending at a 2D-quantum critical point with 
critical exponents z ~ 1 and V ~ 2/3, so the univer- 
sality class of the transition appears to be that of the 
classical 3D-xy model. We have shown that the round- 
ing of the BKT-transition line and the saturation of the 
sheet conductivity close to the QP-transition are remark- 
ably consistent with a gate voltage dependent finite size 
effect. According to this, electrostatic tuning does not 
change the carrier density only but the inhomogeneity 
landscape as well. Taking the resulting finite size effect 
into account we provided consistent evidence for a weakly 
localized insulator separated from the superconducting 
phase by a metallic ground state at quantum critical- 
ity. Consistent with the non-Drude behavior in the nor- 



mal state, characteristics of weak localization have been 
identified in both, the temperature and magnetic field 
dependence of the conductivity. The conductivity along 
the BKT-transition line was found to agree with the stan- 
dard scaling form of quantum critical phenomena, while 
in the weakly localized insulating phase it appears to fail 
in the accessible temperature regime. As in the quan- 
tum scaling approach the scaling function is unknown 
we obtained its form in the superconducting phase. It is 
controlled by the BKT-phase transition line and the vor- 
tex energy. In addition we explored the dependence of 
the vortex core radius and the vortex energy. As the na- 
ture of the metallic ground state at quantum criticality is 
concerned, the limiting dependence of the sheet con- 
ductivity points to Fermi liquid behavior, consistent with 
the evidence for weak localization in the insulating phase 
and non-Drude behavior in the normal state. In conclu- 
sion we have shown that the appearance of metallicity at 
the interface between insulators, a wonderful example of 
how subtle changes in the structure of these systems can 
lead to fundamental changes in physical properties, is a 
source of rich physics in two dimensions. 

This work was partially supported by the Swiss Na- 
tional Science Foundation through the National Center 
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